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ABSTRACT 

We discuss the results from deep Chandra LETG observations of the su- 
pernova remnant 1987A (SNR 1987A). We find that a distribution of shocks, 
spanning the same range of velocities (from ~ 300 to 1700 km s _1 ) as deduced 
in the first part of our analysis (Zhekov et al. 2005, ApJ, 628, L127), can account 
for the entire X-ray spectrum of this object. The post-shock temperature distri- 
bution is bimodal, peaking at kT~ 0.5 and ~ 3 keV. Abundances inferred from 
the X-ray spectrum have values similar to those for the inner circumstellar ring, 
except that the abundances of nitrogen and oxygen are approximately a factor 
of two lower than those inferred from the optical/UV spectrum. The velocity 
of the X-ray emitting plasma has decreased since 1999, apparently because the 
blast wave has entered the main body of the inner circumstellar ring. 

Subject headings: supernova remnants: — supernovae: individual (SNR 1987A) 
- X-rays: ISM 



1. Introduction 

With the rapidly developing impact of the debris of Supernova 1987A with its inner 
circumstellar ring, we have an unprecedented opportunity to witness the birth of a supernova 
remnant, SNR 1987A (McCray 2005). The reappearance of SNR 1987A in radio (Staveley- 
Smith et al. 1992, 1993) and in X-rays (Beuermann, Brandt & Pietsch 1994; Gorenstein, 
Huges & Tucker 1994; Hasinger, Aschenbach & Triimper 1996) at ~ 1200 days after the 
explosion was the first clear sign of this phenomenon. The presence of a triple ring system of 
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relatively dense circumstellar matter (CSM) centered on the supernova (Burrows et al. 1995) 
ensures that astronomers will be able to follow this exciting event for a long time. Indeed, 
we see the continuous brightening in radio (Gaensler et al. 1997; Manchester et al. 2003) 
and X-rays (Park et al. 2004; 2005) as well as the appearance in the optical and ultraviolet 
(UV) of several rapidly brightening spatially unresolved hot spots (Lawrence et al. 2000; 
Michael et al. 2000), which by now encircle the entire inner ring (Sugerman et al. 2002). 

The hydrodynamics of the interaction of the supernova debris with its circumstellar 
matter is quite complex (Chevalier, Blondin & Emmering 1992; Borkowski, Blondin & Mc- 
Cray 1997a,b) and it depends on the details of the density distribution of the CSM as well 
as the density and velocity distributions in the supernova debris. In general, a double shock 
structure forms. The forward shock (blast wave) propagates into the CSM, while a reverse 
shock propagates backwards into the supernova debris. Typically, the separation between 
the blast wave and the reverse shock is ~ 10% of the radius of the blast wave (Chevalier 
1982). Between these two shocks are layers of shocked CSM and shocked supernova debris, 
separated by a (usually unstable) contact discontinuity. These layers are sources of X-ray 
emission. 

Hydrodynamic simulations show that, during the first decades of the supernova remnant 
evolution, the total soft X-ray emission will be dominated by the shocked CSM behind the 
blast wave (Borkowski et al. 1997a). Then, when the blast wave strikes regions of dense gas 
in the inner circumstellar ring, slower shocks will be transmitted into the ring, while reflected 
shocks will propagate backwards through the shocked CSM and will eventually merge with 
the reverse shock. This interaction will further enhance the X-ray emission. By analyzing 
the evolution of the X-ray spectrum of SNR 1987A, we have a unique opportunity to learn 
about the density distribution and elemental abundances of the supernova CSM, as well as 
the physics of high velocity shocks in rarefied gases. These studies may also provide vital 
clues to origin of the inner ring. 

Continuous monitoring of Chandra images of SNR 1987A (e.g. Burrows et al. 2000; 
Park et al. 2002, 2004, 2005; Park et al. 2006, in preparation) have provided us with 
much information about the evolution of the spatial and thermal distribution of the shocked 
gas. However, the poor spectral resolution of the CCD data do not allow us to study the 
kinematics of the X-ray emitting gas in detail. Such information can be obtained only by 
analyzing high-resolution dispersed X-ray spectra. 

The first dispersed X-ray spectrum of SNR 1987A, obtained in October 1999 with the 
Chandra HETG, was discussed by Michael et al. (2002). Unfortunately, the limited photon 
statistics only allowed these authors to construct a composite line profile. The ~ 5000 km 
s _1 FWHM of this profile indicated a blast wave (forward shock) velocity ~ 3500 km s _1 . 
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When we obtained new Chandra (LETG) observations in September 2004, we expected to 
see comparable kinematic velocities; but we were surprised to find that the X-ray emission 
lines had considerably smaller intrinsic widths (~ 500 — 1000 km s _1 ). This observation 
showed clearly that the X-ray emission is now dominated by gas behind shocks of much lower 
velocities than before (Zhekov et al. 2005; hereafter Paper I). 

In Paper I, we measured the widths, line shifts, and fluxes of the strong emission lines 
from hydrogen- and helium-like ions of abundant chemical elements in the spectral energy 
range from 0.4 to 3 keV. (At higher energies the spectrum is dominated by continuum 
emission.) We inferred that the X-ray emission from this object is likely to originate from a 
distribution of shocks with velocities in the range 300 — 1700 km s" 1 . We proposed that this 
ensemble of shocks is a consequence of the blast wave interaction with the complex density 
distribution of the inner circumstellar ring, where transmitted as well as reflected shocks 
play a major role. 

In this, the second part of our analysis of the Chandra LETG observations, we test 
these conclusions further by carrying out a 'global' analysis. In §2 we review briefly our 
observations and data reduction procedures. In §3 we summarize the properties of strong 
X-ray emission lines as described in Paper I. Then, in §4 and §5 we analyze the entire LETG 
spectrum of SNR 1987A and construct a model to fit it. Finally, in §6 we discuss what our 
observations imply about the overall picture of the newborn supernova remnant SNR 1987A 
and its future evolution. 



2. Observations and Data Reduction 

SNR 1987A was observed with Chandra in the configuration LETG-ACIS-S in five con- 
secutive runs (Chandra Obsld: 4640, 4641, 5362, 5363 and 6099) in the period Aug 26 
- Sep 5, 2004 (~ 6398 days after the explosion), providing a total effective exposure of 289 
ksec. The roll angle was chosen so that the negative arm of the dispersion axis was aligned 
approximately north, within ~ 15° from North to West, thus within ~ 10° of the minor axis 
of the inner circumstellar ring (P. A. pa 354° Sugerman et al. 2002). 

We extracted positive (m = +1) and negative (m = — 1) first-order LETG spectra for 
each of the five observations according to the procedure described in the Science Threads for 
Grating Spectroscopy in the CI AO 3.1 1 data analysis software. We merged the resultant 
spectra into one spectrum each for the positive and negative LETG arms with respective total 
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counts of 9,241 and 6,057 in the energy range 0.4 - 7 keV. The difference in photon statistics 
is a result of the different sensitivities of the respective CCD detectors. For each of the five 
data sets we extracted a pulse-height spectrum from the zeroth-order image (according to 
C1AO Science Threads for Image Spectroscopy), contained in a circular region with a radius 
of 10 pixels centered on the source position. We measured the corresponding background 
spectrum from an annulus with inner and outer radii of 10 and 20 pixels, respectively. After 
combining the five pulse- height spectra into one, we measure a total number of 16,557 counts 
in the 0.4 - 7 keV range. We generated the ancillary response functions for all LETG spectra 
(1-st and 0-th order) using the Chandra calibration database CALDB v3.00. 

3. Spectral Lines 

Here we summarize the basic results from our analysis in Paper I of the strong X-ray 
emission lines in the spectrum. We found that the line profiles had no apparent asymmetries 
and that the centroids for the emission lines were consistent with the redshift for the Large 
Magellanic Cloud (see also Fig. 1). Therefore, we could find satisfactory Gaussian fits to the 
line profiles. All the lines seen in (m — —1) LETG spectrum have smaller widths than their 
counterparts in the (m = +1) spectrum. As described in Paper I, this phenomenon results 
from the convolution of the spatial structure and velocity gradients in the X-ray source. 2 It 
is a consequence of the facts that the X-ray source is resolved and that the X-ray emitting 
plasma is confined in the plane of an expanding ring. By having anticipated this physical 
picture and to use this effect to maximum advantage, we chose the roll angle for the LETG 
observations to coincide nearly with the minor axis of the inner circumstellar ring. 

4. Global Spectral Fits 

A systematic approach to modeling the X-ray spectrum of SNR 1987A would be to make 
a global fit of the entire observed spectrum using a model consisting of a distribution of shock 
velocities and ages. Such a global fit has significant advantages over the line ratio analysis 
presented in Paper I: (1) the fit can automatically take into account the quasi-continuum 
due to numerous weak lines; (2) by fitting the shape of the underlying continuum (free-free, 
recombination, and two-photon emission) the model places additional constraints on the 
plasma temperature; (3) the model can constrain the column density of X-ray absorbing 
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gas; and (4) it can yield estimates of relative element abundances. 

It is prohibitive, however, to explore the entire range of models defined by such a 
continuum in parameter space. Instead, we begin by exploring fits to the spectrum with 
a model consisting of one or two plane-parallel shocks and then consider a more complex 
distribution of shocks. 

In adopting this approach, we have used the recent version (11.3.2) of the XSPEC code 
for modeling X-ray emission spectra of shocks. Specifically, we employ the vpshock model, 
which takes into account the non-equilibrium ionization in hot plasmas from version 2.0 of 
NEI models in XSPEC, which are based on ATOMDB (Smith et al. 2001). We augmented 
this atomic database by adding inner-shell processes, which are missing in NEI v2.0. In such 
shock models, the X-ray emission properties are functions of the post-shock temperature 
and the ionization age, n e t, defined as the product of the postshock electron density and the 
time since the gas first entered the shock. Details of these models are found in Borkowski, 
Lyerly, & Reynolds (2001). 

To determine the elemental abundances in our fits, we adopted the same procedure as in 
our previous spectral analyses (Michael et al. 2002; Park et al. 2002, 2004). That is, we only 
varied the abundances of elements having strong emission lines in the observed (0.5-4 keV) 
energy range, namely N, O, Ne, Mg, Si, S and Fe. We fixed the abundances of He and C to 
the values determined by Lundqvist & Fransson (1996) from fits to the optical/UV spectrum 
of the inner ring and the abundances of the remaining elements (Ar, Ca and Ni) to values 
representative of the Large Magellanic Cloud (Russel & Dopita 1992). 

We fitted the positive (m = +1) and negative (m — — 1) LETG 1-st order spectra simul- 
taneously with models having identical plasma characteristics but different line broadening 
parameters defined as in Paper I. Namely, we assumed that the line profiles are Gaussian 
with three sources of line broadening: (1) the spatial extent of the image itself, which can 
be expressed as an equivalent line broadening, independent of wavelength; (2) the thermal 
broadening of the shocked gas; and (3) the broadening due to the bulk motion of the shocked 
gas. For plasma temperatures of interest here (kT as 0.1 — 4 keV) the thermal line broad- 
ening for heavy elements (N through Fe) is negligible compared to the other two sources. 
It is so since the standard relations between the shock parameters yield a ratio of the ion 
thermal velocity to the postshock bulk gas velocity: V th /V bu i k = \f\\ = 0.7 /VA where 
the mean particle weight /i = 0.72 for SNR 1987A and A is the mass of a given ion in units 
of the proton mass. Therefore, the thermal velocity of various ionic species is only between 
10% and 20% of the bulk gas velocity (A = 14 [N] through 56 [Fe]). On the other hand, the 
spatial size of SNR 1987A in wavelength units (AA = 0.047 A, Paper I) results in a line 
broadening larger than that from the bulk gas velocity of the shocked gas. Therefore, we 
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express the net line width (FWHM) as: 

A\ tot = 2AA ± 2z (A/A ) Q A , (1) 

where the plus (minus) sign refers to the m — +1 (m — — 1) spectrum, respectively. The 
first term on the right hand side of equation (1) represents the broadening due to the spatial 
extent of the source. The second term represents the broadening due to the bulk motion. 
The power-law function of wavelength in the second term allows for the possibility that the 
mean bulk velocity of shocked gas emitting a given line may depend on the excitation or 
ionization stage of the emitting ion. The parameter zq determines the line broadening at 
some fiducial wavelength, A , and the power-law index a is to be determined. 



4.1. Discrete Shock Models 

We fitted one- and two-shock models to the LETG spectra of SNR 1987A (rebinned to 
have a minimum of 30 counts per bin). The results from our model fits are shown in Table 1. 
We find that a fit with one-shock model is statistically unacceptable and that it yields an 
unreasonably small X-ray absorption column density. (This value is from 3 to 4 times smaller 
than that derived from the previous X-ray observations, Michael et al. 2002; Park et al. 2004; 
whose results were consistent with the neutral hydrogen column density towards SNR 1987A 
deduced from observations in the optical and UV, Fitzpatrick & Walborn 1990; Scuderi et 
al. 1996). Figure 4 also illustrates how difficult is for the one-shock model to give a good 
fit to sectral lines of ionic species which require different plasma conditions (e.g., Si XIII 
and Si XIV with ionization potential, IP, of 2.67 and 2.44 keV and Fe XVII with IP of 1.27 
keV). Therefore, we favor the two-shock model. The results for the line broadening from the 
global fit are close to those derived in Paper I: z = 1.13 x 10 -3 [7.4 x 10~ 4 — 1.54 x 10~ 3 ] vs. 
1.57 x 10- 3 ; A = 0.061 [0.057 - 0.064] vs. 0.047 A; and a = -1.17 [-1.94 - -0.08] vs. 
-1.3 (the 90%-confidence limits are given in brackets). The relatively small differences can 
be attributed to the fact that the global fitting procedure simultaneously handles both lines 
and continuum. 

As a check, we also fitted one- and two-shock models to the pulse-height spectrum of the 
zero-order LETG image. The best-fit (two-shock) results are shown in Table 1. Of course, the 
parameters determined from the pulse height spectrum, especially the abundances, are less 
tightly constrained than those determined from the dispersed images. However, agreement of 
the fit to the pulse-height spectrum with the fit to the dispersed spectrum gives us confidence 
in using also the pulse-height spectrum to infer element abundances and conditions in the 
X-ray emitting gas. 
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The abundances of SNR 1987A derived here are consistent with those derived from 
previous CCD spectra (Michael et al. 2002, Park et al. 2004). We find that the abundances 
of Ne, Mg, Si, S and Fe are very close to those found in X-rays for the LMC SNRs (Hughes 
et al. 1998). On the other hand, the N and O abundances derived here are lower by a factor 
of about two than those found for the inner ring by Lundqvist & Fransson (1996). We note 
that similar lower nitrogen abundance was found by Pun et al. (2002) in their analysis of the 
optical and UV spectrum of Spot 1 on the inner circumstellar ring. But, it is interesting to 
note that the relative abundance of nitrogen and oxigen derived in our analysis, N/O = 1.1 
(by number), is consistent with that derived from the modeling of the optical/UV emission 
of the inner ring (Lundqvist & Fransson 1996, Sonneborn et al. 1997). Thus, we confirm 
the nitrogen enhancement in the CSM around SNR 1987A, which strongly suggests that 
the CSM consists of matter that had undergone CNO processing in progenitor star and was 
subsequently ejected. 

Significant departures from electron-ion temperature equilibration may be present in the 
post-shock plasma (Michael et al 2002). In the fits we present here, we have not included 
this effect. We have explored models including departures from electron-ion temperature 
equilibration (the vnpshock option in XSPEC) and we find that they give equally good fits 
to the spectra, but a detailed analysis of such models is beyond the scope of this paper (see 
also §6). 



4.2. Distribution of Shocks 

One of the basic results from the line-profile analysis of Paper I was that a distribution 
of shocks with velocities in the range 340 — 1700 km s _1 is compatible with the LETG data. 
Given the relation between the postshock temperature and the shock velocity for strong 
adiabatic shocks with electron-ion temperature equilibration, 

kT e = = 1.4[y s /1000 km s" 1 ] 2 keV; (2) 

16 

where \i = 0.72m p for SNR 1987A. The observed velocity range corresponds to a temperature 
range of kT e = 0.15 — 4 keV. 

To develop this scenario, we constructed a new XSPEC model having the following fea- 
tures: (1) the distribution of emission measures in the shocked gas is determined from the 
Chebyshev polynomial algorithm as is used in the standard XSPEC cQpvmkl model (Lemen 
et al. 1989); (2) the basis vectors for X-ray emission from the distribution of shocks are 
those from the XSPEC vpshock model; (3) all shocks share the same element abundances 
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and X-ray absorption. Moreover, to decrease the number of free parameters, we also as- 
sumed: (4) the ionization age (n e t)i of each individual shock % is a power-law function of the 
postshock temperature, i.e., {n e t)i oc (kT e )^ where p is a free parameter. Details about the 
line broadening are given below. 

As in our earlier analysis, we assume the emission lines to be Gaussian profiles with 
FWHM depending on the wavelength. We recall that the wavelength dependence of the line 
width represents the physical notion that faster shocks produce higher-temperature plasma 
whose emission is dominated by lines at shorter wavelengths. Now, by fitting the entire 
spectrum with models having a distribution of shock velocities, we can test this assumption 
in a more rigorous way. Namely, we assume that the total X-ray spectrum of a parcel of 
shocked plasma (having temperature T e ) is moving with postshock velocity given by equation 
(2). Then, similarly to eq. (1) we write: 

A\ tot = 2AA ± 2z s \ , (3) 

where z s = (3(3/A)V s /c ; V s is the shock velocity; (3/4) stands for the case of a strong shock 
entering a stationary gas; and c is the speed of light. The parameter f3 (0 < f3 < 1) allows 
for various geometrical effects such as viewing angle to be considered as well. 

With these assumptions, we can find a satisfactory (x 2 /dof = 414/431) simultaneous fit 
to the positive and negative LETG spectra of SNR 1987A. Figure 2 shows the distribution 
function of emission measures and ionization ages determined by the fit. Note that the shape 
of the inferred shock distribution is bimodal peaking very close to the temperature values 
derived in our discrete two-shock model fit. To verify this conclusion, we tried to fit various 
models having more X-ray plasma with temperature between the two peaks. In no case 
could we find an acceptable fit without a bimodal distribution. Note also that the ionization 
age of the shocked gas decreases with increasing shock temperature. 

The absorption column density and element abundances have the values (given are la 
errors): N H = 1.52 ± 0.02 x 10 21 cm- 2 , N= 0.88 ± 0.13, 0= 0.10 ± 0.01, Ne= 0.31 ± 0.02, 
Mg= 0.25 ± 0.02, Si= 0.29 ± 0.02, S= 0.48 ± 0.09, and Fe= 0.16 ± 0.02. We note that 
these parameters are also very close to those determined from the two-shock model fit (for 
comparison see Table 1). Moreover, the derived value of f3 — 0.35 ± 0.06 is close to that 
expected for a radial gas motion in a disk viewed at inclination angle i: sin sin i = 0.45 
(sin0 = 2/tc , < <fi < 7r/2; i = 45° for the inner ring; Sugerman et al. 2002). These results 
give us confidence in our global fits to the LETG spectrum of SNR 1987A. 
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5. MARX simulations 

As a final test of our model, we used the MARX 3 software to simulate the actual 
two-dimensional images and dispersed spectra which incorporate the kinematic and spectral 
information deduced from our analysis of the line profiles, line ratios and global fits to the 
observed X-ray spectra of SNR 1987A. Since the object is spatially resolved by Chandra, 
our simulations take into account the actual brightness distribution of the X-ray emission as 
well. The simulations correspond to the actual observing conditions (e.g. date of observation, 
exposure time, pointing of the telescope, roll angle etc.). 

The actual brightness distribution of the X-ray emission used in our simulations is shown 
in Fig. 3. This is the 0-th order image from our LETG observations deconvolved according 
to the same procedures used in our imaging observations (Burrows et al. 2000; Park et al. 
2002, 2004). Shown also are the eight sectors we used to represent the details of the image 
morphology. 

To construct the input spectra, we assume that the X-ray emitting plasma is confined in 
the plane of the inner ring and that it is expanding radially. Locally, we assume a distribution 
of shocks with characteristics as described in §4.2 (see also Fig. 2), with emission measure 
scaled to correspond to the brightness distribution in each sector. 

We constructed our MARX simulations according to the following procedure. For 
each sector, we simulated the observed spectrum, assuming that the X-ray spectrum corre- 
sponds to the physical picture given above and a flux scaled to the total observed flux from 
SNR 1987A according to the fractional brightness of that sector in the zero order image. We 
determined the red- or blue-shift of the spectrum of each sector according to its postshock 
plasma velocity, corrected for its position (azimuthal angle), and the inclination of the inner 
ring (i = 45°). We then concatenated the resultant simulated data from the all sectors to 
derive a data set for the entire object. We then extracted the simulated LETG spectra 
following the CIAO Science Threads for Grating Spectroscopy. In order to minimize the 
statistical fluctuations in the simulated spectra, we averaged the ten different realizations of 
the simulated spectra. 

Figures 5, 6 and 7 display these average spectra overlaid with the observed positive and 
negative first order LETG spectra. The excellent correspondence between the simulated and 
actual spectra gives us confidence that our model is a reasonable approximation to the actual 
interaction of the blast wave with the inner ring. 
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6. Discussion and Conclusions 

The analysis of the LETG spectra of SNR 1987A presented in §4 and §5 of this paper 
confirms and refines the main conclusions of Paper I. The X-ray emitting plasma is largely 
confined to a radially expanding ring near the inner optical ring. The X-ray emitting gas 
is heated by a distribution of shocks and is not in ionization equilibrium. The electron 
temperature distribution of the shocked gas is bimodal, with peaks at ~ 0.5 and ~ 3 keV, 
respectively. 

Our analysis leads us to a general picture of SNR 1987A in which the X-ray emission 
is dominated by two components. The hot (kT e ~ 3 keV) component results from shocks 
of relatively high velocity that propagate through relatively low density (n e ~ 10 2 cm -3 ) 
circumstellar gas. The cool (kT e ~ 0.5 keV) component results from shocks transmitted into 
the denser (n e ~ 10 4 cm -3 ) gas of the inner circumstellar ring. 

The fact that the inferred ionization age decreases with shock velocity (Fig. 2) is con- 
sistent with this picture. The slower shocks propagate through the higher density gas, while 
the faster shocks propagate through the lower density gas. Therefore, if the typical timescale 
since the shock has entered the gas is comparable (t ~ few years), the ionization age, n e t, 
of the gas behind the slower shocks will be substantially greater than that of the gas behind 
the faster shocks. 

In their analysis of a spectrum of SNR 1987A taken with the High Energy Transmission 
Grating of Chandra in October 1999, Michael et al. (2002) inferred an electron temperature, 
kT e ~ 2.6 keV from the line ratios. By stacking all the emission lines into a composite 
profile, they also a measured the width, FWHM ps 5000 km s _1 , of a composite line profile, 
from which they inferred a 'typical' shock velocity ~ 3500 km s" 1 that was consistent with 
the radial expansion velocity of the X-ray image as measured by Park et al. (2002). 

Since the gas behind such a shock would have electron temperature ~ 17 keV if 
the electrons and ions were in temperature equilibrium, Michael et al. concluded that 
the electron and ion temperatures had not equilibrated. This was a reasonable scenario, 
since the electron-ion temperature equilibration time, (n e t) = 4 x 10 u (kT keV ) 1 ' 5 cm" 3 s 
(Spitzer 1962), of such a gas would be substantially greater than the typical ionization age, 
n e t ~4x 10 10 cm -3 s, of shocked gas having n e 400 cm -3 and t ~ 3 years. 

But the shock environment of SNR 1987A has changed considerably since that observa- 
tion. The line profiles of the current (September 2004) observations (Paper I) indicate shock 
velocities 300 — 1700 km s" 1 , and hence typical post-shock temperatures kT e ~ 0.1 — 4 keV. 
Moreover, such low shock velocities suggest that the shocked gas must have greater density, 
say, 4 x 10 4 > n e > 1600 cm" 3 . From these considerations and Fig. 2, we conclude that 
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there is no longer any compelling evidence that the electrons and ions in the shocked gas are 
out of temperature equilibrium, although this may still be the case for gas behind the faster 
shocks propagating through the lower density gas. 

We expect the rapid evolution of the X-ray spectrum to continue. Park et al. (2005) 
found that the X-ray light curve is brightening at an accelerating rate. This fact and the 
correlation of the X-ray image with the optical hotspots (Park et al. 2006, in preparation) 
indicate that the blast wave is interacting with much denser gas than in the past. This 
observation is supported by the noticeable deceleration of the radial expansion velocity of 
the X-ray image. The most recent Chandra data show that this velocity has decreased from 
~ 4000 km s^ 1 (Park et al. 2004) to ~ 1600 km s^ 1 (Racusin et al. 2006, in preparation). 

The rapid brightening of the X-ray light curve of SNR 1987A (Park et al. 2005) ensures 
that we may soon be able to obtain spectra with excellent photon statistics. Thus, with 
future Chandra observations we should be able to constrain the model further and subject 
it to some critical tests. 

When the blast wave strikes dense clumps of gas in the circumstellar ring, the interaction 
will give rise to both transmitted and reflected shocks, both of which will cause enhanced 
X-ray emission correlated with the optical hotspots. In both cases the velocity distribution 
of the shocked gas will shift toward lower velocities. Therefore, we expect that the X-ray 
emission lines will continue to become narrower. 

The evolution of the temperature distribution of the X-ray emitting gas will tell us the 
relative contribution of reflected and transmitted shocks to the rapidly brightening source. 
If the gas behind the transmitted shocks dominates the X-ray emission, the bimodal tem- 
perature distribution of the X-ray emitting gas will become more skewed towards lower 
temperatures. In fact, Park et al. (2005) have already noticed this trend in their two-shock 
model fits to Chandra nondispersed (ACIS-S) spectra over the past five years. On the other 
hand, if the gas behind the reflected shocks contributes substantially to the X-ray emission, 
the hot component in the X-ray spectrum will continue to remain substantial. 

The radio emission from SNR 1987A is non-thermal (Gaensler et al. 1997; Manchester et 
al. 2002) and is probably produced by relativistic electrons accelerated by the reverse shock 
(Manchester et al. 2005). Park et al. (2005) have noted that the hard X-ray light curve 
resembles the radio light curve of SNR 1987A and have suggested that the radio emission 
and the hard X-rays might have a common origin. This suggestion raises the question of 
whether the X-ray spectrum may have a non-thermal (i.e., power-law) component. Up to 
now, we have seen no compelling evidence for the presence of such a component, either in 
the pulse-height spectra (Michael et al. 2002) or in the LETG spectrum discussed here. On 
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the other hand, we cannot rule out the possibility of a significant non-thermal component, 
especially at photon energies > 4 keV. 

Observations of broad line profiles of L a and H a (Michael et al. 2003; Smith et al. 
2005) show that the supernova debris is crossing this reverse shock with velocities ~ 10 4 km 
s _1 . If the hard X-rays originate from the region behind the reverse shock and are produced 
by thermal processes, then we would expect to see strong emission lines of high-ionization 
species with relatively broad profiles (~ few 10 3 km s _1 ). On the other hand, if the hard 
X-rays are predominantly non-thermal, we would expect to see a decrease in the equivalent 
widths of the X-ray emission lines at higher energy. 

Future Chandra observartions will be also helpful in resolving the issue with the CNO 
abundances derived in X-rays. As discussed in § 4 (see also Table 1), the N/O ratio is 
consistent with the value derived from the analysis of the optical/UV spectra of the inner 
ring in SNR 1987A, but the total amount of light metals (C, N, O) derived in our analysis 
(CNO pa 0.15 solar) is a factor of two lower than that derived from optical/UV spectra (CNO 
w 0.28 solar, see Lundqvist & Fransson 1996). The latter ratio is consistent with the average 
CNO abundances for LMC (CNO ps 0.26 solar, Russel & Dopita 1992). This means that the 
X-ray analysis confirms the enhancement of the N/O ratio, which may be attributed to CNO 
processing of the matter by the progenitor star. But the X-ray observations indicate a deficit 
in the net abundance of C + N + O (the CNO processing cannot alter their total amount). 
On the other hand, the abundances of heavier metals (Ne, Mg, Si, S and Fe) derived here 
are consistent with typical values for the SNRs in LMC (Hughes et al. 1998). 

What could account for the low abundances of C + N + O derived from the X-ray 
spectrum? On the one hand, we cannot rule out the possibility of purely technical reasons 
for this discrepancy, such as greater uncertainties in the X-ray spectrum at low energies 
due to the CCD degradation or uncertainties in the atomic data. On the other hand, there 
might be physical reasons. Could the amount of CNO elements in the progenitor star be 
considerably lower (a factor of two) than the average for LMC? Or, could there be an extra 
source of (likely nonthermal) continuum, the effect of which would be to reduce the CNO 
abundances if interpreted as thermal continuum? Future Chandra grating observations with 
anticipated better photon statistics will help us find the correct answer. In addition, such 
observations may also reveal some evolution (increase) for the abundances of elements like 
O, Mg, Si and Fe resulting from the destruction of dust from the supernovae debris and the 
inner circumstellar ring by the shocked hot gas emitting in X-rays. 

This work was supported by NASA through Chandra Awards G04-5072A (to CU, Boul- 
der, CO) and GO4-5072B (to NCSU, Raleigh, NC). 
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Table 1: Shock Fitting Results (1-st and 0-th order spectra) 



Parameter 




LETG 1-st 




LETG 0-th 




1-shock 




2-shock 




2-shock 




X 2 /dof 


661/436 




417/433 




121/111 




Ntf(10 21 cm- 2 ) 


0.63 [0.45 - 


0.821 

J 


1.47 [1.25 - 


1.73] 


2.13 [1.52 - 


• 3.00] 


kTi (keV) 






0.50 [0.47 - 


0.54] 


0.41 [ 0.35 


- 0.66] 


kT 2 (keV) 


2.16 [2.08 - 


2.21] 


2.72 [2.44 - 


3.02] 


3.40 [2.63 ■ 


- 4.061 


taui^ 






5.21 [3.73 - 


7.83] 


6.36 [3.10 - 


- 29.11 


tau 2 ( a ) 


0.69 [0.62 - 


0.77] 


1.43 [0.77 - 


15.4] 


1.43 [0.76 - 


• 13.51 


EM/ 6 ) 






5.51 




12.7 




EM 2 ( fe ) 


2.45 




1.69 




1.44 




H 


1 




1 




1 




He(2.57) 


2.57 




2.57 




2.57 




C (0.09) 


0.09 




0.09 




0.09 




N (1.63) 


0.40 [0.32 - 


0.50] 


0.77 [0.60 - 


1.00] 


0.12 [0.0 - 


0.73] 


(0.18) 


0.055 [0.047 - 0.066] 


0.092 [0.075 - 0.104] 


0.10 [0.073 


- 0.15] 


Ne(0.29) 


0.25 [0.22 - 


0.28] 


0.29 [0.25 - 


0.34] 


0.29 [0.21 ■ 


- 0.41] 


Mg(0.32) 


0.21 [0.18 - 


0.24] 


0.24 [0.20 - 


0.28] 


0.23 [0.16 ■ 


• 0.33] 


Si(0.31) 


0.24 [0.21 - 


0.27] 


0.28 [0.22 - 


0.32] 


0.66 [0.47 ■ 


- 1.12] 


S (0.36) 


0.42 [0.28 - 


0.56] 


0.45 [0.32 - 


0.60] 


0.43 [0.24 ■ 


- 0.96] 


Ar(0.54) 


0.54 




0.54 




0.54 




Ca(0.34) 


0.34 




0.34 




0.34 




Fe(0.22) 


0.16 [0.14 - 


0.18] 


0.16 [0.14 - 


0.17] 


0.11 [0.08 - 


- 0.18] 


Ni(0.62) 


0.62 




0.62 




0.62 




(0.5-2 keV) 


1.42 




1.47 




1.55 




F^ (0.5-6 keV) 


1.80 




1.84 




1.95 





The 90%-confidcnce intervals are given in brackets and all abundances are expressed as ratios to their solar 
values (Anders & Grevesse 1989). For comparison, the inner-ring abundances of He, C, N, and O (Lundqvist 
& Fransson 1996); those of Nc, Mg, Si, S, and Fe typical for the LMC SNRs (Hughes et al. 1998); and 
Ar, Ca and Ni abundances representative for LMC (Russell & Dopita 1992) are given in the first column in 
parentheses. Note that the H, He, C, Ar, Ca, and Ni abundances were kept fixed in all model fits (see text). 
a (n e t) in units of 10 11 s cm -3 . 

b EM — J n e riHdV in units of 10 58 cm' 3 , assuming a distance of 50 kpc. 
c X-ray flux in units of 1CP 12 ergs cm -2 s _1 . 
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Fig. 1. — Left panel: Line shifts for the strong emission lines in the positive (diamonds) 
and negative (squares) LETG first-order spectra of SN 1987A. The solid line in the left panel 
represents the average red-shift for the Large Magellanic Cloud (z = 9.53 x 10~ 4 ). Right 
panel: The OVIII L a profile in the positive (solid line) and negative (dashed line) LETG 
first-order spectra. The count rate in the m — — 1 spectrum has been multiplied by a factor 
1.5 to correct for the lower sensitivity of the detector in that arm. 
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0.1 1.0 

Temperature [keV] 

Fig. 2. — Emission measure (EM) distribution of shocks in SNR 1987A as derived from the 
DS model with 25 points logarithmically spaced in the (0.1-4 keV) postshock temperature 
range. The two vertical dashed lines indicate the plasma temperature derived from the 
discrete two-shock model. The solid line shows the derived ionization age of each shock 
(n e t). 



Fig. 3. — The deconvolved 0-th order image from the LETG observations of SNR 1987A. 
North is up and East is on the left. Shown are the boundaries of the eight sectors (subimages) 
used in the MARX simulations (§5). The horizontal line (1 arcsec) indicates the scale of the 
image. 
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12.5 13.0 13.5 14.0 



13.0 13.5 




4.5 15.0 15.5 16.0 16.5 17.0 17.5 




4.5 15.0 15.5 16.0 16.5 17.0 17.5 





19.0 19.2 19.4 



19.0 19.2 19.4 



Fig. 4. — The positive 1-st order LETG X-ray spectrum of SNR 1987A (rebinned to have a 
minimum of 30 cts per bin) near some strong spectral lines : (i) empty squares with la-error 
bars - observed spectrum; (ii) solid line - fits from the XSPEC 1-shock (left column) and 
2-shock (right column) models (see Table 1). Horizontal axes - observed wavelength (A); 
vertical axes - flux density (photons s^ 1 A -1 ). 
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Si XIV (+1 ) 




Si XIV (-1) 




6.00 6.10 6.20 6.30 





6.40 6.50 6.60 6.70 6.80 6.90 7.00 



6.40 6.50 6.60 6.70 6.80 6.90 7.00 



Mg XII (+ 1 ) 




0.006 r 
0.005 r 

0.004 - 

0.003 r 

0.002 
0.001 f- 
0.000 



Mg XII (- 1 ) 




.20 8.30 8.40 8.50 8.60 




Mg XI (- 1 ) 




9.00 9.10 9.20 9.30 9.40 



Fig. 5. — The X-ray spectrum of SNR 1987A near the strong spectral lines of the H-like 
and He-like ions Si XIV, Si XIII and Mg XII, Mg XI: (i) empty squares with lcr-error bars - 
observed spectrum; (ii) solid line - simulated MARX spectrum. Horizontal axes - observed 
wavelength (A); vertical axes - flux density (photons s _1 A -1 ). The positive and negative 
1-st order LETG spectra are shown in the left and right columns, respectively. 
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Nc X (+1) 




Nc X (-1 ) 



0.01 4 r- 

0.01 2 7 
0.01 7 
0.008 r 
0.006 1 



0.004 - , 

: 1 rT- 

0.002 - 



0.000 




2.00 12.10 12.20 12.30 



12.00 12.10 12.20 12.30 




13.2 13.4 13.6 13.8 14.0 




13.2 13.4 



Fc XVII (+1 ) 





0.000 

14.80 14.90 15.00 15.10 15.20 15.30 15.40 



Fc XVII (-1 ) 




0.000 

14.80 14.90 15.00 15.10 15.20 15.30 15.40 





Fig. 6. — The same as in Fig. 5 but for the H-like and He-like ions iof Ne X and Ne IX and 
the strong lines of Fe XVII and VIII (Ly^). 
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Fc XVII (-1 ) 



ft, 



16.6 16 



7.0 17.2 17.4 



16.6 16.8 17. ( 



17.2 17.4 



VIII (+1 ) 




0.006 p- 

0.005 r 

0.004 - 

0.003 r 
0.002 r 



O VIII (-1 ) 



9.00 19.10 19.20 



18.80 18.90 19.00 19.10 19.20 




O VII (-1 ) 



21.4 21.6 21.8 22.0 22.2 22.4 



21.4 21.6 21.8 22.0 22.2 22.4 




N VII (-1 ) 



24.50 24.60 24.70 24.80 24.90 25.00 25.10 



24.50 24.60 24.70 24.80 24.90 25.00 25.10 



Fig. 7. — The same as in Fig. 5 but for Fe XVII, the H-like and He-like ions iof O VIII (Ly a ) 
and O VII and the H-like ion of N VII. 



